Study on Evolvement Complexity in an Artificial Stock Market 
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An artificial stock market is established based on multi-agent . Each agent has a limit memory of 
the history of stock price, and will choose an action according to his memory and trading strategy. 
The trading strategy of each agent evolves ceaselessly as a result of self-teaching mechanism. Simu- 
lation results exhibit that large events are frequent in the fluctuation of the stock price generated by 
the present model when compared with a normal process, and the price returns distribution is Levy 
distribution in the central part followed by an approximately exponential truncation. In addition, 
by defining a variable to gauge the "evolvement complexity" of this system, we have found a phase 
cross-over from simple-phase to complex-phase along with the increase of the number of individuals, 
which may be a ubiquitous phenomenon in multifarious real-life systems. 

PACS numbers: 89.90.-|-n,02.50.Le,64.60.Cn,87.10.-|-e 
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In social systems, such as insect societies, increased 
colony size is associated with profound and wide-ranging 
changes in "internal" organization and operation 1, '3|. 
For instance, larger colony size is correlated with in- 
creased homeostasis, cooperative activity, spatial orga- 
nization of work, and caste polymorphism to name but 
a few "social correlates" [l|, Q • Gautrais et al catch up 
with a model demostrating a proximate mechanisms to 
emerge polymorphism in insect societies, which indicates 
that specialization only occurs above a critical colony 
size such that smaller colonies contain a set of undiffer- 
entiated equally inactive individuals while larger colonies 
contain both active specialists and inactive generalists, as 
has been found in empirical studiesjSj. The specialization 
of workers upon certain tasks can increase colony pro- 
ductivity. The experimentation Weidenmuller made in- 
dicated that the dynamics of the colony response changed 
as colony size increased: colonies responded faster to per- 
turbations of their environment when they were large (60 
or more individuals) than when they were small|j|. These 
findings provide intriguing new examples of the ways in 
which individuals, each using only local information, act- 
ing simply and independently and not subject to any cen- 
tral or hierarchical control, can coordinate group-level be- 
havior which differs from that of each individual, as does 
economical system. Every economical agent behaves sim- 
ply contrast to the system which is composed of them. 
Economical system complex behaviors also result from 
repeated nonlinear interaction between each others. But, 
does the social economical system have similarity to in- 
sect societies that macro-properties have to do with par- 
ticipants size? In this letter, we have found a phase cross- 
over from simple-phase to complex-phase along with the 
increase of the number of individuals based on our model, 
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which may be a ubiquitous phenomenon in multifarious 
real-life systems. 

There are many modelling methods to explain origins 
of the observed behavior of market price as emerging 
from simple behavioral rules of a large number of het- 
erogeneous market participants, such as behavior- mind 
modeip, [a, d ynamic-games model Q , multi-agent model 
H [lalll 111 111 El and so on. The mainstream 
method is agent-based modelling because of its simple- 
ness, agility and verisimilitude and which based on a styl- 
ized description for the behavior of agents. Here, we pro- 
posed a stock market model based on multi-agent that 
incorporates the feedback between the price trend and 
agent's trading strategy. Therefore, our model will de- 
mostrate that each agent has a limit memory of the his- 
tory of stock price and will choose an action according to 
his memory, and that the trading strategy of each agent 
evolving ceaselessly as a result of self-teaching mechanism 
will influence the price trend inversely, which resemble 
the minority game|l5J. 

In our model, before a trade, each agent should choose 
an action: to buy, to sell or to ride the fence, the for- 
mer two should determine the price and amount of the 
trading-application. The buyer with higher price and the 
seller with lower price will trade preferentially, and the 
trading-price is the average of selling-price and buying- 
price. The stock price is the weighted average of trading- 
price according to the corresponding trading-amount 16J. 

Each agent holds a so-called decision-matrix, which 
can tell him how to do according to the history of stock 
price. Let p{t) be the stock price at time step i, then the 
range of fluctuation is 

/(t) = {p[t) - p{t - l))/p{t - 1) e (-1, +oo) (1) 

For the sake of simplification, the range of fluctua- 
tion is categorized into 5 types: drastic fall (/ € 
(-1,-0.05)), faU (/ e [-0.05,-0.01]), near immovabil- 
ity( / G (-0.01,0.01)), rise (/ e [0.01,0.05]) and drastic 



TABLE I: John's decision matrix. A. sell all his shares; B. 
sell half of his shares; C. do nothing; D. spend all his cash on 
shares; E. spend half his cash on shares 

Patterns Action A Action B Action C Action D Action E 
0.40 0.20 0.10 0.20 



rise (/ £ (0.05, +oo)), which are denoted by -2, -1, 0, 1 
and 2 respectively. The agent's memory is limited to the 
current 5 fluctuations, thus there are 5^ = 3125 different 
fluctuation-patterns. The decision-matrix contains the 
probabihties of trading strategies according to the dif- 
ferent fluctuation-patterns. For instance, table 1 shows 
a decision-matrix of an agent named John. Based on 
this matrix, John will choose to sell half of his shares in 
hand at probability 0.40 when the present fluctuation- 
pattern is "1, 0, -1, 0, 2". If an agent decides to buy 
or to sell, the buying-price or selling-price will be cho- 
sen completely randomly in the interval [p{t), l.lp(i)] or 
[0.9p{t),p{t)] respectively. 

After a trade, each decision-matrix will change as a 
result of self-teaching mechanism. For each agent, if his 
action made his money increase, the corresponding prob- 
ability in his decision-matrix will be doubled, contrarily, 
it will be halved. After that, the probabilities under the 
very pattern will be normalized. For instance, if John's 
action were an unsuccessful one, the probabilities under 
the fluctuation-pattern "1, 0, -1, 0, 2" would become 
"0.125, 0.25, 0.25, 0.125, 0.25" after normalization. Ap- 
parently, there will be no changes if the agent did noth- 
ing or his action kept his money unaltered. In order to 
mimic the "bounded rationality" and "inductive think- 
ing" of investors (Ili^ 18] , we set a very small probability 
7, which is called the reversal parameter. Agents may 
change their decision-matrix in completely contrary di- 
rection at the probability 7. 

When proper initial condition and parameters have 
been chosen, the artificial stock market can generate its 
stock price. In figure 1, we report a typical simulation 
result about price time series generated by our model. 
Which is similar to the reality (inset). In this simula- 
tion we set the market size as 1000 (i.e. 1000 stockhold- 
ers), the initial stock price as 50. The initial quantity of 
fund and shares owned follows uniform distribution in the 
interval [0, 1000000] and [0, 10000] respectively, and the 
original fluctuation-pattern are randomly selected from 
the 3125 candidates. Notice that an agent's action may 
be restricted by his wealth. In other words, he may be 
prevented from buying or selling because of, respectively, 
a shortage of fond or shares in hand. 

Large numbers of simulations have been performed to 
check if the model can generate price time series of key 
characteristics according with the reality. Since chaotic 
characteristic is one of complex dynamical properties of 
economical system, which has been demonstrated by pre- 
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FIG. 1: Time series of the typical evolution of the stock price 
in the interval t € [0, 40000], where 7 — 0.01 and the elements 
of initial decision-matrices are chosen completely randomly 
in the range [0, 1] before normalization. The insect is the 
Dow Jones Industrial Average(DJIA) from 01-02-1931 to 12- 
31-1987. 



vious studies[T9l l2(l]. We have calculated the Lyapunov 
exponent and correlative dimension of the stock price 
time series generated by the model, carried out princi- 
pal component analysis, and drawn the conclusion that 
our model can not only create stock price trends rather 
similar to the real, but also show the chaotic behavior in 
deep consistency with the real stock market. The details 
are omitted, and can be referred to the corresponding 
reference |2l|. 

In addition, we have calculated the distribution of price 
return r{t), where r(i) is deflned as the difference between 
two successive logarithms of the price: 



r{t) = p{logp{t + At) - logp{t)) 



(2) 



Here, /3 is a positive constant. The corresponding price 
returns with At — 2 are shown in flgure 2, from which 
one can see that large events are frequent in the fluctu- 
ation of the stock price generated by the present model 
when compared with a normal process, which agrees with 
the previous empirical studies well 22, 23, 24, 25, 26]. In 
figure 3, we plot the probability distribution of price re- 
turns, and the fitted Gaussian curve for the case Ai = 1. 
Comparing with normal distribution, the present returns 
distribution is of more peaked center and fatter tail, ac- 
cording with the empirical studies that suggest the dis- 
tribution of returns in real- life financial market is a Levy 
distribution in the central part followed by an approxi- 
mately exponential truncation 22, 23, 24, 25, 26]. Since 
our main goal in this letter is not to show the comparison 
between price time series generated by our model and the 
reality, more details are omitted here. 

In succession, let's discuss whether macro-properties 
have to do with participants size. As far as whether the 
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FIG. 2: The price returns with At = 2, where /3 = 10^ and 
< t < 10000. 



FIG. 4: 
10000, ; 
[0,0.1]. 



The evolvement complexity of size A^, where L = 
I — 100 and 7 is randomly chosen in the interval 
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FIG. 3: The normalized probability density of price returns 
over different time scale At = 1, 2, 4, 8, 16, 32. The solid curve 
is the fitted Guassean curve for the case At = 1. 



system is complex and how m.uch is the complexity degree 
are concerned, there isn't a criterion. Whether the sys- 
tem is complex is sometimes apparent. This, however, is 
not necessarily the case. If there is only one simple time 
series generated by a certain system, it is very hard to 
make sure that the system is complex. On the contrary, 
we can easily find the evidence if the system is obvious 
incomplex. In order to answer that problem , we loosely 
define a variable, i.e. evolvement complexity under 
the condition that the behaviors of our artificial econom- 
ical system are similar to those of the reality. In gen- 
eral, our model will run long time. But the model with 
small market size will run more shorter time than that 
with large market size and will generate price time series 
which tail is monotone or belongs to a finite period or 



reaches a approximate fixed point. Such a price time se- 
ries p{t) generated by our model with length L is defined 
to be simple. Here, "tail" means the last O.IL points, and 
"it reaches an approximate fixed point" means the ratio 
of its range to its average is smaller than 0.05. Let the 
market size be fixed and other parameters be variable, 
if m simple time series are generated by n independent 
experiments, then the evolvement complexity of size N 
is loosely defined as C{N) = 2^^. Note that the length 
of tail (O.IL) and the criterion of what is an approximate 
fixed point (< 0.05) are not especially selected to gen- 
erate the simulation results shown in the present letter. 
One can write a program and easily check that the simu- 
lation results are robustious for a wide parameter space, 
thus the following phenomena are. 

As is shown in figure 4, the system evolutive behavior 
is evidently divided into 3 areas. When N < 50, the be- 
havior is simple, and its evolvement complexity increases 
rather slowly with the increasing of the number of indi- 
viduals. When N > 100, the system has a great evolve- 
ment complexity, but its "complex degree" does not in- 
crease with the increasing of the number of individuals. 
Therefore the areas TV < 50 and A'' > 100 can be consid- 
ered as the simple-phase and complex-phase respectively. 
Between these, the complexity increases fiercely with the 
increasing of the number of individuals, thus it is called 
the critical interval with the inf-critical-point 50 and sup- 
critical-point 100. Here, we have found a phenomenon 
that the behavior of our artificial economical system is 
correlated with increased participants size, which is sim- 
ilar to that of insect societies. 

The economical system constitutes one among many 
other systems exhibiting a complex organization and dy- 
namics with similar behavior, which, with large num- 
ber of mutually interacting parts, self-organize their dy- 
namics with novel and sometimes surprising macroscopic 



("emergent") properties. The phenomenon that the 
evolvement complexity can be divided into 3 parts ow- 
ing to the increasing of the number of individuals maybe 
one of the common characteristics among various com- 
plex systems that do not seem alike at all in appearance. 
The conception "evolvement complexity" is novel and in- 
teresting, but it is hard for us to give out a strict and 
appropriate definition. Therefore, it is a innovative point 



as well as a shortcoming in this letter. Although the def- 
inition is rough, we believe it will enlighten physicists on 
how to measure the complexity of complex systems. 
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